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THE EFFECT OF THE PARAMETERS OF THE MOLDING UNIT 
IN GLASS-MOLDING MACHINES ON THE NONSTATIONARY 
TEMPERATURE FIELD OF GLASS MOLDS 

V. I. Volkov 1 

Translated from Steklo i Keramika, No. 12, pp. 11 - 15, December, 1998. 

The temperature field of glass molds in the case of their external cooling is considered. The temperature field 
is represented as a sum of two components: the stationary component not depending on time, and the compo¬ 
nent cyclically varying with time. The cyclic component is investigated for the case where the heat flows at 
the mold boundaries are harmonic functions of time. Ratios establishing the relationship between the mold 
temperature variations per cycle and the molding unit parameters, as well as the relationships between the op¬ 
timum values of these parameters, are obtained. The amplitude of temperature variations on the mold working 
surface is calculated for different parameters of the molding unit. 


Thermocyclic loads in various fields of engineering con¬ 
stitute a basic factor which determines or has a perceptible 
effect on equipment wear. Glass and casting molds, cylinders 
in internal combustion engines, and cutting tools can serve as 
examples [1]. Therefore, the solution of the problem of in¬ 
creasing the service life of this equipment is directly related 
to the problem of determination of the nonstationary temper¬ 
ature fields arising inside this equipment in the course of op¬ 
eration. 

The solution of this problem is especially important for 
the glass industry, since fluctuations in temperature on the 
working surface of the mold not only determine the longev¬ 
ity of the mold, but the quality of molded articles as well and 
the very possibility of glass molding. The less the tempera¬ 
ture fluctuations, the higher the quality of the glass product, 
since in this case the molding proceeds under more homoge¬ 
neous conditions. Moreover, a decrease in the amplitude of 
temperature fluctuations makes it possible to raise the aver¬ 
age temperature on the molding surface (without the danger 
of reaching the temperature of glass adhesion), which also 
improves the quality of glass articles [4, 5]. The amplitude of 
the mold temperature fluctuations is a factor limiting the effi¬ 
ciency of glass-molding machines, since with an increase in 
the thermal load on a mold, the temperature fluctuations on 
its working surface inevitably increase, and if their value ex¬ 
ceeds the difference between the glass adhesion point and the 
temperature producing notches, “hammering,” and other de¬ 
fects, the process of glass molding becomes impossible. 


1 Orient Ltd., Moscow. 


In order to determine the nonstationary temperature field 
of the glass mold, we decompose temperature t( r, r) into two 
components: the stationary one u (r) which does not depend 
on time, and the periodically varying with time v(r, t): 

/(r,T) = M(r) + v(r,x), (1) 

where r are space coordinates and x is time, and consider 
each of them [6]. This approach is based on the following 
simple transformation: Jet Q( t) be the heat flow density, and 
Q* the average heat flow density on the mold working sur¬ 
face per cycle (assuming that the heat flow on this surface is 
homogeneous, i.e., does not depend on its coordinates). Ob¬ 
viously, 

1 r 

Q*=j- | 0(T)rfT, 

1 0 

where T is the molding cycle duration, i.e., the time interval 
between the penetration of two consecutive glass melt drops 
(for a parison mold) or two consecutive bullets (for a finish¬ 
ing mold) into the mold. 

Let us denote the difference Q (x) - Q* by the function 
q (x); then we see that the heat flow on the working surface 
of the mold is decomposed into two components: the station¬ 
ary Q* and the periodically changing with time q (x), when 

0(x)=2* + <7(x). (2) 
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The described transformations are illustrated in Fig. 1 for 
the step function Q(x): 


0(x) = 


Q v ai nT<*x< nT + BT 

0,at BT + nT <x<(n + \)T,n =0,1,2,... 


where n is the number of the period; Q t is the heat flow den¬ 
sity; B is part of the cycle during which the glass mixture 
stays in the mold. 

Similarly, for the outer surface of the mold: represent the 
thermal efficiency coefficient a(x) as the sum 

a(x) = ct* + a 0 (x), (3) 

where a* is the average coefficient of thermal efficiency on 

the mold outer surface per cycle, when 
T 

a* = \/T Ja(x)</x, and a Q (x) is equal to the difference 
o 

a(t) - a*. 

Substituting Eq. (1) into the heat conduction equation, 
and Eqs. (2) and (3) into the respective boundary conditions, 
we obtain two systems of equations, one of which, for u (r ), 
does not contain the variable x at all. The stationary compo¬ 
nent of the mold temperature n(r ) is well studied [6]; there¬ 
fore we will specify only the equations for determining the 
variable component v(r, x). This is the heat conduction equa¬ 
tion 


dv 

—=aAv, 
dx 


( 4 ) 


where a is the temperature conductance of the mold material; 
A is the Laplace operator, and the boundary conditions 


dr • 


( 5 ) 


dv 

-X-l-^ix + x^uir^) 


1 + - 




+ a * v. 


( 6 ) 




a 
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Fig. 1. Heat flow density on the mold working surface Q decom¬ 
posed into the stationary Q* and the cyclic q components: a) density 
of heat flow Q(x) and the average heat flow density per cycle Q* 
(stationary component); b ) cyclic component of heat flow density. 


waves do nor affect each other due to the high degree of ab¬ 
sorption. If the mold wall is sufficiently thin, the amplitude 
of the resulting wave, owing to the interference, can become 
either higher or lower than the initial wave amplitudes. 

The above model does not take into account heat transfer 
at the ends of the mold, and, consequently, the higher the 
amount of thermal energy abstracted from the molded glass 
article which is diffused on the mold outer surface, the more 
precisely will the heat transfer process in the mold be de¬ 
scribed. As was shown both in the experiments [6] and in in¬ 
dustrial tests [7, 8], only a relatively small (~ 20%) part of 
the energy is abstracted through thermal radiation and heat 
transfer in the sites of the mold fastening. Therefore, the heat 
flow inside the mold wall can perceptibly differ from the ra¬ 
dial heat flow only within a limited area near the upper and 
bottom ends of the mold. 

Consider the simplest case where the variable compo¬ 
nents of the heat flows at the mold boundaries are described 
by harmonic functions depending on time. Let, for example. 


where X is the heat conductance of the mold; t a is the average 
temperature of the cooling agent (air) surrounding the outer 
surface of the mold [6]. 

The time interval t 0 takes into account the possibility of 
the periodic processes of mold heating and cooling being 
asynchronous. For example, at x 0 >0, the cooling of the 
mold begins x 0 seconds prior to the moment in which glass 
melt (for the parison mold) or the bullet (for the finishing 
mold) gets into the mold. Thus, the variable heat flows at the 
mold boundaries excite the heat waves which propagate to¬ 
ward each other inside the mold wall. If the wall is thick, the 


q(x) = q x sin cox; 

(x 0 (x + x 0 ) = a, sin (cox + e) , 

where co = 2n/T is frequency; e = cox 0 is the phase difference 
between the processes of mold cooling and heating. 

This case, which is the simplest for heat flows varying 
periodically with time, is at the same time a fundamental 
one, since any periodic functions q (x) and a 0 (x + x 0 ) can be 
decomposed into the Fourier series by sinus and cosines 
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Fig. 2. Dependence of the optimum cooling phase e 0 on the mold 
wall thickness l. 


(provided they are integrable and at least piecewise continu¬ 
ous, which is always fulfilled for the actual heat flows). 

Equation (6) can be simplified by ignoring the term 
v/(u-( a )<l. The measurements performed on an IS-6 - 2 
machine producing half-liter bottles of dark green glass [9] 
showed that the ratio between the temperature variations am¬ 
plitude V and the mean temperature u on the outer surface of 
the mold is ~ 6%. 

By approximating the cylindrical wall of the glass mold 
by a plane parallel wall, the solutions of Eqs. (4) - (6) can be 
expressed by elementary functions. Assuming that x 0 = 0, 
and x, = / (/ is the thickness of the mold wall), we find the 
variable temperature component for a depth of x, v(x, t): 

v (x, t) = V(x ) sin [tax + arg(F(x ))], (7) 

where arg (V (x)) and V are the argument and modulus, re¬ 
spectively, of the complex amplitude V equal to: 


Fig. 3. Variable component of temperature on the mold working 
surface: 1 — A ~ 1, e = e 0 ; 2 — A = 0.71, e = e 0 ; 3 — A = 0; 4 — 
A = 0.71, £ = e 0 + tt; .5 — A = 1, e = e 0 + 7t. 

molds cooled by the verti-flow system). Therefore, in the 
computation of Eq. (7), the terms containing Bi* can be ig¬ 
nored. 

The amplitude of temperature variations on the mold 
working surface V 0 is of practical interest, since its signifi¬ 
cant increase can bring about the production of defective ar¬ 
ticles and acceleration of the mold wear due to increased 
thermocyclic loads. As seen from Eq. (8), the value of the 
amplitude V depends on the mold wall thickness / and the 
cooling phase e. Denote by e„ the cooling phase in which V 0 
reaches the minimum level for a given value of /. By ignor¬ 
ing the members containing Bi*, we obtain from Eq. (8) 

e 0 = arctan (tan kl th kl ). (9) 


ch [(1+ /)£(/-x)] + Bi*sh[(l + /)£(/-x)3~/l exp(/s) 

* sh [(1 +m + Bi * ch [(1+ i)kf\ * ( 

where i is the imaginary unit, i 2 = - 1; = V<o/2 a is the wave 

CL * 

number; Bi* =- is the modified Bio number; A = 

(1 + i)Xk 

a i [«(/)-*,] 

9i 

In our case. A - 1, since the heat flows on the mold sur¬ 
faces are equal. 

The real values of | Bi* | in glass molding are not high. 
Thus, for the 2W-12-01 machine equipped with a cooling 
system ensuring the most intense heat removal from the 
mold compared to other domestic glass-molding machines, 
a* « 500 W/(m 2 • K) [6], and assuming T= 10 sec, we get 
| Bi | » 0.05. For other domestic machines | Bi* | will be 
even lower, and for foreign machines with a closed cooling 
system, this value will be similar (we are not discussing here 


The plot of the function e 0 (/ ) is shown in Fig. 2. Fig¬ 
ure 3 shows the dependence of the amplitude of temperature 
variations on the working surface on the mold wall thickness 
for e = e 0 . As can be seen, with a sufficiently thick wall, the 
amplitude V 0 starting with a certain value of /, (when 
kl t a 5.5) does not depend at all on the mold wall thickness. 
Indeed, the temperature waves excited on the mold surfaces 
do not affect each other, the mold wall can be regarded as 
semi-infinite, and we arrive at the known [10] relationship 

v= q l /(-j2kX)e~ k * sin % j • 

At / < /,, the interference of the temperature waves 
causes a monotonic decrease in the temperature variation 
amplitude on the molding surface, as the wall thickness in¬ 
creases (provided e = e 0 ), up to V 0 = 0 at / = 0 (A = 1). If the 
cooling phase e differs from e 0 , the amplitude V 0 with a de¬ 
crease in / will decrease more slowly than at e = e 0 , and its 
minimum will be nonzero and will be attained at a certain 
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Fig. 4. Variable component of the mold working surface tempera¬ 
ture: / — A = 3.65, £ = e 0 ; 2 — A = 1.41, e = £ 0 ; 3 — A = 0; 4 — 
A= 1.41, e = £ 0 + tc; 5 — A = 3.65, £ = £„ + it. 

value / = / 0 * 0. A further decrease in l <l 0 causes a rapid in¬ 
crease in V 0 at / —> 0. 

If e significantly differs from c 0 , V 0 can monotonically 
grow with the decrease of /. In this case (for example, when 
cooling is delayed), before the “cooling wave” excited at the 
outer surface of the mold reaches the molding surface, the 
latter will be heated to a significantly high temperature (as 
/ 0), close to the temperature of molded glass. This heat¬ 

ing will be more intense the thinner the mold wall (and ac¬ 
cordingly, the lower the mass and the heat capacity of the 
mold). The “cooling wave” which arrives virtually without 
absorption due to the low value of the wall thickness will sta¬ 
bilize the process, and yet, due to the delay, the cooling will 
last longer than required, which will cause a significant de¬ 
crease in the temperature of the working surface (as l-> 0) to 
a level close to the initial temperature of the cooling agent. 
Obviously, the thinner the mold wall, the more intense is the 
cooling process (Fig. 3). 

In the above described model of the heat transfer taking 
place in the mold in the course of glass molding, the parame¬ 
ter/! was equal to 1. However, in the more general case, the 
heat processes at the mold boundaries can be described not 
by harmonic but by more complex periodic functions of 
time. Then, as these functions are decomposed into a Fou¬ 
rier’s series, the amplitudes of the harmonics of one order at 
the opposite surfaces of the mold will not be equal and, ac¬ 
cordingly, the parameter^ can take an arbitrary value. In this 
case, the shape of the curve describing the dependence of the 



Fig, 5. Dependence between the parameter A and the wall thickness 
I at which the temperature wave amplitude minimum on the mold 
working surface is attained: 1 ) A (/ 0 ) at A < \;2)A(l 0 )atA> 1. 


amplitude of the working surface temperature variation V 0 on 
the mold wall thickness will be changed. 

With A < 1 (A being a positive value), same as with 
A = 1, the amplitude of the variable component of the mold 
temperature, starting with a certain value of /, depending on 
A , decreases, as / decreases, at e = e 0 . However, the curves 
have a minimum distinct from zero at / 0 * 0, and the value of 
this minimum and the value of / 0 are both functions of the pa¬ 
rameter A. Figure3 shows the dependences V 0 (l) for 
A = 0.71 (curves 2 and 4 ) and A = 0 (curve 3 ). The latter 
means that the mold is constantly chilled, i.e., a(x) = const. 
Note that the specific shape of the curves V(l) is determined 
by the parameter A: the greater is A, the wider and deeper the 
curve minimum, and the higher the value of /, starting with 
which the amplitude of temperature variations on the mold 
surface decreases to a minimum at £ = e 0 . And vice versa, at 
£ = e 0 + 7i, the greater A, the sooner the temperature varia¬ 
tions increase as the mold wall thickness increases. For all 
other values of £ the curves V(l ) lie within a region which is 
limited on one side by the “minimum” curve corresponding 
to e = e 0 and on the other side by the “maximum” curve cor¬ 
responding to £ = e 0 + n, for example, between curves 1 and 
5 in Fig. 3 at A = 1, and between curves 2 and 4 for A = 0.71 

If A > 1, the “minimum” curves F 0 (l ) will always have a 
minimum equal to zero at a certain value of /„ depending on 
A, i.e., K 0 (/ 0 ) = 0. At the same time, while for A, < 1 the 
value of / 0 decreases with an increase in A up to / 0 = 0 at 
A = 1, for A > 1, on the contrary, /„ increases with an increase 
in A, and curve F 0 (/) has the specific shape at the minimum 
point (Fig. 4). 

The plot establishing a relationship between the parame¬ 
ters A and / 0 is shown in Fig. 5. For a sufficiently large A (in 
fact, for A > 3) the value of / 0 can be approximately found as 
l^\/(k\n2A). 

For the purpose of our calculation, the glass mold wall 
was regarded as plane-parallel. Let us define the area of ap¬ 
plication of this approximation. It is well known that a cylin¬ 
drical wall in thermal calculations can be approximated by a 
plane-parallel plate if the cylindric surface radius is signifi¬ 
cantly larger than the wall thickness, i.e., r Q » 1. Since the re- 
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gion of effective decrease in the temperature variation ampli¬ 
tude corresponds to values of the mold thickness below 1 cm 
(Fig. 3 and Eq. (8)), the above cited relationships are true at 
least for the finishing molds used for bottle production, and 
even more so for the molds used to produce glass jars. 

It seems reasonable to offer another criterion for the as¬ 
sessment of suitability of the above model for description of 
temperature variations inside a cylindrical wall. Taking into 
account that a decrease or increase in the temperature varia¬ 
tion amplitude depending on the value of c is determined by 
the effect of the temperature wave interference, the better the 
cylindrical temperature waves propagating in an actual mold 
can be approximated by the plane waves, the more precisely 
the results obtained for a flat wall will describe the processes 
taking place in the mold. Such approximation will be justi¬ 
fied if the width of the wave-front significantly exceeds its 
length, i.e., A c 2nr fl (A is the length of the temperature 
wave, 2jtr 0 is the width of the temperature wave-front). For a 
cast iron mold at T= 10 sec, A ~ 1.2 cm. Therefore, the crite¬ 
rion of applicability of the plane-parallel approximation will 
be r 0 /Vr >0.5 x 10' 3 for r 0 expressed in meters and T ex¬ 
pressed in seconds. 

The proposed approach to the description of glass mold 
temperature makes it possible to determine uniquely such 
significant molding parameters as the beginning of the mold 
cooling process (for periodic mold cooling) and the thickness 
of the mold wall. The obtained relationships can be accepted 
as a first approximation to precise results, since the heat 
flows at the mold boundaries were assumed to be harmonic 
functions of time. In order to use them in engineering de¬ 
signs, the effect of the type of time dependence of the speci¬ 
fied heat flows on the temperature field should be investi¬ 


gated. However, being periodic functions of time, the heat 
flows at the mold boundaries can be decomposed into a Fou¬ 
rier series, i.e., represented as sums of harmonic functions. 
Thus, the obtained results can be used as instruments in such 
investigation. 
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